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MAGNETIC DIFFUSIVITY TENSOR AND DYNAMO EFFECTS IN ROTATING AND SHEARING TURBULENCE 

A. Brandenburg', K.-H. Radler-, M. Rheinhardt^, and R J. Kapyla'-' 
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ABSTRACT 

The turbulent magnetic diffusivity tensor is determined in the presence of rotation or shear. The question is 
addressed whether dynamo action from the shear-current effect can explain large-scale magnetic field genera- 
tion found in simulations with shear For this purpose a set of evolution equations for the response to imposed 
test fields is solved with turbulent and mean motions calculated from the momentum and continuity equations. 
The corresponding results for the electromotive force are used to calculate turbulent transport coefficients. The 
diagonal components of the turbulent magnetic diffusivity tensor are found to be very close together, but their 
values increase slightly with increasing shear and decrease with increasing rotation rate. In the presence of 
shear, the sign of the two off-diagonal components of the turbulent magnetic diffusion tensor is the same and 
opposite to the sign of the shear. This implies that dynamo action from the shear-current effect is impossible, 
except perhaps for high magnetic Reynolds numbers. However, even though there is no alpha effect on the av- 
erage, the components of the a tensor display Gaussian fluctuations around zero. These fluctuations are strong 
enough to drive an incoherent alpha-shear dynamo. The incoherent shear-current effect, on the other hand, is 
found to be subdominant. 
Subject headings: MHD - turbulence 



1. INTRODUCTION 

Many of the stellar and planetary magnetic fields are be- 
lieved to be the result of a dynamo process that converts ki- 
netic energy from turbulent motions and shear into magnetic 
energy. A particular challenge consists in explaining the field 
on length scales that exceed the scale of the turbulence. This 
topic has traditionally been addressed within the framework 
of mean-field electrodynamics (Krause & Radler 1980). 

Over the decades the applicability of this theory has re- 
peatedly been questioned (e.g., Piddington 1981, Vainshtein 
& Cattaneo 1992). Meanwhile, direct simulations of hy- 
dromagnetic turbulence have begun to show dynamo action 
(Meneguzzi et al. 1981, Meneguzzi & Pouquet 1989, Nord- 
lund et al. 1992, Brandenburg et al. 1996, Cattaneo 1999). 
In some particular cases, large-scale fields are being gener- 
ated (Glatzmaier & Roberts 1995, Brandenburg et al. 1995, 
Brandenburg 2001) which raises the question about the mech- 
anism responsible for this phenomenon. In cases where the 
flow is systematically non-mirror symmetric the association 
with an a effect is obvious. However, there are now also 
examples of nonhelical large-scale dynamos owing to turbu- 
lence under the influence of shear alone (Brandenburg 2005a, 
Yousef et al. 2007). Their interpretation is not straightfor- 
ward, because several possible mechanisms have been pro- 
posed that might produce dynamo action from turbulence and 
shear alone, i.e. without rotation and stratification that other- 
wise would have been the main ingredients of an a effect. The 
most detailed investigations have been carried out in connec- 
tion with the so-called shear-current effect (Rogachevskii & 
Kleeorin 2003, 2004, Radler & Stepanov 2006, Rudiger & 
Kitchatinov 2006). Another possibility is a magnetic a ef- 
fect that is driven by a current helicity flux, as was suggested 
by Vishniac & Cho (2001; see also Brandenburg & Subrama- 
nian 2005c). A third possibility might be an incoherent (ran- 
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dom) a effect with zero mean and finite variance, suggested 
by Vishniac & Brandenburg (1997) in connection with accre- 
tion discs (see also Sokolov 1997, Silant'ev 2000, Fedotov et 
al. 2006, Proctor 2007). The only reliable way to determine 
what is the dominant effect is to calculate all relevant compo- 
nents of the a and turbulent magnetic diffusivity tensors in a 
general expansion of the electromotive force in terms of the 
mean magnetic field. 

The case considered in Brandenburg (2005a) is unnecessar- 
ily complicated because the shear employed there depends on 
two Cartesian coordinates. A simpler possibility is to consider 
a shear flow depending linearly on only one coordinate and we 
shall pursue this idea in the present paper The shear-current 
effect and the incoherent a effect could then still operate. Be- 
cause we will use periodic boundary conditions there can be 
no magnetic helicity flux, so the Vishniac & Cho (2001) ef- 
fect is then ruled out, even though it could still, at least in 
principle, explain the generation of a mean magnetic field in 
the simulations of Brandenburg (2005a), which do possess a 
helicity flux. 

In this paper we calculate all relevant components of a,j 
and Tjijic using the so-called test field method. This method 
was introduced by Schrinner et al. (2005, 2007) in connection 
with convection in a spherical shell and used later by Branden- 
burg (2005b), Sur et al. (2007) and Brandenburg et al. (2008) 
in connection with forced turbulence in Cartesian boxes. The 
essence of this method consists in solving evolution equations 
for the fluctuations of the magnetic field around suitably de- 
fined test fields such that all relevant coefficients can be com- 
puted. 

2. BASIC ASSUMPTIONS AND METHOD 

In the following we introduce first the mean electromotive 
force and its relation to the mean magnetic field. We then dis- 
cuss the equations describing the turbulent flow that eventu- 
ally leads to this electromotive force and explain the test field 
method used to calculate the coefficients which relate it to the 
mean field. Particular attention is paid to the possibility that 
the shear-current effect may lead to self-excitation of mean 
magnetic fields. 
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2.1. The turbulent electromotive force 
In mean-field electrodynamics the behavior of the mean 
magnetic field B depends crucially on the mean electromo- 
tive force £> - uxb, where u and b denote the deviations of 
the fluid velocity U and the magnetic field B from their mean 
parts U and B, respectively. For sufficiently weak variations 
of B in space and time, and if there is no small-scale dynamo 
producing a mean electromotive force on its own, we have 



&, = aijBj + TjijkdBj/dxk 



(1) 



with tensors aij and rjijk determined by u and U. 
In this section it is sufficient to define mean quantities like 

£ or B, referring to Cartesian coordinates (x, y, z), simply by 
averaging over all x and y. Below, in ^2.31 a different defini- 
tion will be introduced that covers and refines this simple one. 
Clearly B can now no longer depend on x and y and hence all 
its first-order derivatives can be expressed by the components 
of V X B = i-dBy/dz, dBJdz, 0). SUghtly deviating from the 
usual notation, in which V x B is equal to fioJ where fio is the 
vacuum permeability, we put in this paper simply V x B = 7, 
being aware that / is then no longer exactly the electric cur- 
rent density. Instead of ([T]i we may now write 



£; = OijBj - JJijJj 



(2) 



with a new tensor 77,7 defined such that 77,1 - t/qs and 77,2 = 
-77,13. As 73 = the 77,3 are without interest and we may put 
them equal to zero. 

We further consider the background turbulence, which oc- 
curs in the absence of rotation or shear, as homogeneous, 
isotropic and mirror-symmetric. Then we have even under 
the influence of rotation or shear Oij - 0; see Appendix lAl 

As for T]ij consider first the case of rotation of the fluid with 
an angular velocity Q, which defines the Coriolis and cen- 
trifugal forces and is assumed to be aligned with the z-axis. 

The actual turbulence is then again homogeneous but no 
longer isotropic. Instead it is axisymmetric with respect to the 
z-axis, that is, all mean quantities depending on the turbulent 
velocity field are invariant under arbitrary rotations about the 
z-axis. We may then conclude by usual symmetry arguments 
that 

= m^ij + SeijkClk + 6'CliClj , (3) 

where Q - Q/Q with Q = |£2|, and 770, 6 as well as 6' are 
spatially constant coefficients, which may depend on D.. So 
we arrive at _ _ _ 

£ = -770/ + 6D.XJ. (4) 

Since 73 = the 6' term in ([3]) is without influence. The last 
term in equation (HJi describes the ilxj effect (Radler 1969). 
Whereas 770 approaches a nonzero value as Q ^ (the value 
determined by the background turbulence), 6 vanishes like Q. 
Note that = 0. 
Consider next the case with shear defined by the velocity 
= (Q,S X, 0). Now the actual turbulence is again homoge- 
neous but no longer axisymmetric. In view of the application 
of symmetry arguments we consider first in the more gen- 
eral (coordinate-independent) form S g(h ■ x) where g and h 
are unit vectors which are orthogonal to each other, and x is 
the position vector. The only available construction elements 
for rjij are then 6ki, etim, g and h, for due to the homogeneity 
of the turbulence, 77,j cannot depend on x. Thus we have 

Tlij = m^U + Kigigj + Kihihj + K-igihj + K4gjhi + ■■■ (5) 



where 770, k\, K2, kj, and are spatially constant coefficients 
and the dots stand for additional terms containing etim- The 
aforementioned coefficients may depend on 5. (A depen- 
dence on scalars defined by g and h is without interest since 
g^ - - \ and g ■ h - 0.) The terms containing ej^/,„ have 
structures like eijkgk or giCjkigkhi. Since S g{h ■ x) is invari- 
ant under simultaneous sign changes of g and h, rjij must also 
have this property and so these terms have to be cancelled. 
Returning now to - {0,Sx,0), that is g = (0, 1,0) and 
/i = (1,0,0), we see that 



with 



£ - -r]oJ - K - J 



(Kn K12 0\ 
K21 K22 




(6) 



(7) 



This covers the "shear-current effect" (Rogachevskii & Klee- 
orin 2003). We may assume that 770 is independent of S (that 
is, it is determined by the background turbulence alone). Then 
K\i and K22 are even functions of S that vanish like 5^ as 
5 — > 0, whereas k\2 and K2\ are odd functions that vanish 
like S . Again we have £- = 0. 

In both cases, with rotation or with shear, we may restrict 
our attention to 



< 2 



(8) 



£i = -^ijJj' ^^iJ 

The four quantities 77,^ are simply related to 770 and 6, or 770 
and the k/j respectively. 

2.2. Turbulence with rotation or shear 

We consider a compressible fluid satisfying an isothermal 
equation of state. In the absence of rotation or shear the mo- 
mentum and continuity equations can be written in the form 



dU 

dt 



— = -[/■ Vl/-<Vlnp-H/-HFv 



d\np 
dt 



-f/- Vlnp-V-f/, 



(9) 



(10) 



where c, is the sound speed, here considered as constant, p the 
mass density and / a random forcing function. Furthermore, 
^'visc = p^'V ■ 2pvS is the viscous force, and S,j - \{Uij + 
Ujj) - ■ U is the traceless rate of strain tensor. 

To come as close as possible to the assumptions on the 
background turbulence adopted above, i.e., homogeneity, 
isotropy and mirror-symmetry, the forcing function / was 
specified for a cubic domain of size L x L x L as follows. 
During each time-step / is a single transverse (solenoidal) 
plane wave proportional to x e where the wavevector is 
taken randomly from a set of pre-defined vectors with compo- 
nents being integer multiples of In/L and moduli in a certain 
interval around an average value which we simply denote by 
kf, and e is an arbitrary random unit vector not aligned with 
^f. The corresponding scale, In/kf, is also referred to as the 
energy-carrying scale of the turbulence. Moreover, the time 
dependence of / is designed to mimic ^-correlation, which 
is a simple and commonly used form of random driving (cf. 
Brandenburg 2001). Nevertheless, owing to inertia, the cor- 
relation time of the turbulent velocity is of course finite, even 
for perfect (5-correlation of the forcing. 

As mentioned above, we simulate the turbulence in a fi- 
nite domain using (shearing-) periodic boundary conditions 
02. 51 ). Then the background turbulence can be at no instant in 



a strict sense homogeneous, isotropic or mirror-symmetric. It 
would approach these properties if the ratio of the size of this 
domain and the scale of the forcing function (that is fcf/fci with 
^1 = InjL) became very large. There are, however, practical 
bounds on this ratio. For moderate values, which we have to 
accept, the background turbulence approaches the mentioned 
properties only after averaging over long times. Then, of 
course, the turbulence appears also as statistically steady. By 
these reasons mean quantities, that is, averages over x and y, 
which are derived from the turbulence, show still fluctuations 
in z and f, and these disappear after averaging over sufficiently 
long time intervals. 

When rotation is added, two new terms arise on the right 
hand side of equation (|9]l, the Coriolis force, -2£i x U, and 
the centrifugal force, (Q x jc) x Q. The latter is unimportant 
for weak compressibility, to which we restrict ourselves in the 
following, and this term would also not be compatible with 
periodic boundary conditions, so it is neglected. 

Turning now to the case with shear we redefine the velocity 
U by splitting off the shear term , that is {7 -> t/ + [/^ . This 
implies 



vu + u vu^ + 



VU\ (11) 



The second term on the right hand side corresponds to an ad- 
ditional advection with the mean flow and will be subsumed 
in the definition of an advective derivative. 



D/Dt = d/dt + Sxd/dy. 



(12) 



The third term is equal to S U^y, where y is the unit vector in 
the y direction. The last term in equation ( fTTT l vanishes. Thus, 
equations (|9]l and (fTOl i turn into 



^ ^-U-VU + SUJ-c^lnp + f + F, 



and 



D\np 
Dt 



= -U -Vlnp-V -U. 



(13) 



(14) 



It should be noted that U resulting from these equations 
is not purely turbulent, but also contains a large-scale flow 
which provides an additional shear and therefore a mean vor- 
ticity. This is qualitatively suggestive of a hydrodynamic 
mean-field effect analogous to the shear-current effect; see 
Elperin et al. (2003). 

In this paper we deal, apart from one exception, with the 
purely kinematic problem, so there is no Lorentz force in 
equation ( fT3] ). In 34. 11 the fully nonlinear problem is consid- 
ered and hence the Lorentz force is included in the momentum 
equation. 

2.3. Test field method 

Proceeding now to consequences of the induction equation 
we consider primarily the case of shear, in which the fluid 
velocity is U + . In the case of rotation we have to put 
equal to zero. We further represent B according to B = V x A 
by a vector potential A. Uncurling the induction equation and 
using a suitable gauge transformation of A we find'* 



DA 

— = -SAyX + UxB-T]J . 



(15) 



t/j(Aj,-A,j). 

The second term is an advection term and the first term can be written as 
U^Aji = -Uj-Aj plus a gradient term that can be removed with a gauge 

transformation. Note also that U^j-Aj = SAyXj. 



This equation as well as those derived from it in what follows 
apply to the case of rotation if D/Dt is replaced by d/dt and 
S is put equal to zero. 
Now we define a mean field F belonging to the field F as 

F{x,y,z)^-j F{x + ^,y + ri,z)d^dT] (16) 

J-L/2 J-L/2 

with L as specified above. The following comments on the 
definition in equation ( fTSI l as well as equations (fTTl i and ( fTSl l 
below apply, however, even if L is an arbitrary length, not nec- 
essarily related to the domain size. Our definition ( fT6l ) implies 
that averaging of F commutes with taking any derivatives of 
F with respect to x, y, z or f, that is, the sequence of these op- 
erations can be changed. In what follows we also use the rule 

FG - F G, which applies exactly if F is independent of x and 
y, and has otherwise to be considered as an approximation. It 
is however only needed in cases in which that independence 
of X and y can be justified, that is, in which it applies exactly. 
Clearly, F is independent of x and y if F is periodic in x and y 
with the period length L. We note further that, owing to (fT6b . 
we have H - x, and that therefore has to be considered as 
a mean field. 
Taking now the average of ( fTSl l we obtain 

^jy^ _ _ 

— = -SAyX + UxB + uxb-j]J. (17) 



In view of the determinationof fi = « x ft we are interested in 
b - V X a, where a - A - A. Taking the difference between 
equations (fTSl l and ( [TtI i we obtain 

Da 



— — = -Sa,.x + Uxb + uxB + uxb-uxb-rii, (18) 
Dt ■ iJ^ y J 

where j - J - J. 
In order to determine the quantities rjij introduced above we 

specify B in the relevant relations such that it is equal to one 
of the elements out of a set of test fields, B'', and denote the 
corresponding £, J , etc., by Efl, respectively. Then, in 
particular, equation ^ turns into 



(19) 



After having calculated the numerically for two properly 
chosen B"^ we may then determine the four rjij. 
For the calculation of the we apply ( fTsT i ^, 

_ _ 

— = -S4x + Vxb'' + UXB'' +uxb'' -uxb''-T]f . (20) 

Although an additional mean flow can develop in some of the 
simulations with shear (see above), this term is still weak and 
is neglected in the following, hence we put 17 = in equa- 
tion ( |20] |. 

As test fields B'' we may use, e.g., the fields B**^ defined by 
B'" = B(cosfe,0,0), B-" = B(0,cosfe,0) (21) 

with a constant B and a constant wavenumber k. Denoting the 
corresponding £fl by £fl'^ we find 

= rii2Bk sin kz , £■ ' = -///i Bk sin kz, i^l,2. (22) 

^ In the corresponding eq. (27) of Brandenburg (2005b), the U term is 
incorrect. However this did not affect his results because U either vanished 
or it consisted only of a shearing motion that was treated correctly in the code. 
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After having calculated the fi**^ these equations allow us to 
determine the rjij. In order to avoid difficulties at the zeros of 
sin kz it is useful to carry out the calculations with test fields 
defined analogously to the B'''^ but with sinfe instead 
of cos kz. For the corresponding £^J"* we find then equations 
analogous to ( |22] | but with - cos kz instead of sin kz. From 
this and equation (l22l i we obtain immediately 

i]a = -{Bk)-^{&f sin kz - Sf cos kz) 

T]i2 = ' sin kz - cos kz), / = 1 , 2 (23) 

(see also Brandenburg 2005b). 

We recall that for homogeneous turbulence, which is con- 
sidered here, the t/,- ,- ha ve to be independent of z. That is, the 
cos kz and sin kz in i22l l and ( |23] l should be compensated by z- 
dependencies of the £f!'' and However, due to fluctuations 
(cf. ^2.2l i. no perfect compensation can be expected. 

We further recall that in ([U and so also in (|2]i all derivatives 
of B that are higher than first-order have been ignored. By 
this reason the results for the rjij obtained with the above test 
fields apply exactly only in the limit A: — > 0. In general there 
is a dependence of the rjjj on k. This corresponds to a non- 
local connection between £> and B, which is considered here 
only in a very weak sense (by taking into account first-order 
derivatives of B). In a more general sense it is investigated in 
Brandenburg et al. (2008). Here we have used k = ki where 
^1 means the smallest finite wavenumber in the z-direction in 
the domain in which the turbulence is simulated, ki - In/L; 
see 

The results for 77,^ are also independent of the value of B. 
If one wanted to address the question of nonlinearity, which 
is not the purpose of this paper, one must also solve equa- 
tion (fTsT i and allow the resulting magnetic field to feed back 
onto the flow via the Lorentz force. 

For the discussion of the results concerning 77,^ we introduce 
the quantities 

rjt ^ j(rin + mi) , + e ^ \{m\ - mi) ■ (24) 

In the case of rotation we put further 

d^^imi-mx) (25) 

and expect e to be equal to zero, while 7712 and 7721 have to 
have the same nonzero moduli, but opposite signs so that 5 
is nonzero. With shear, however, e can be in general nonzero 
and there is no simple relation between 7/12 and 7721. 

It is convenient to present tji in normalized form and express 
it in terms of the quantity 



?7to 



1 "rms'^f 



(26) 



3 "rms'^f ' 

which corresponds to the result for 771 obtained under the first 
order smoothing approximation applied to the high conduc- 
tivity limit under the assumption that the correlation time is 
given by {u^mf^kf)^^, i.e., that the Strouhal number is unity (cf. 
Brandenburg & Subramanian 2005b). 

2.4. Dispersion relation 

In the case of rotation without shear there are only decay- 
ing solutions of the mean-field equations. This can be easily 
seen from the energy balance equation for the mean magnetc 
field (see Radler 1980). The situation with shear alone is how- 
ever different and the possibility of a so-called shear-current 
dynamo is still under debate (see Rogachevskii & Kleeorin 



2003, 2004, Brandenburg 2005b, Radler & Stepanov 2006, 
Riidiger & Kitchatinov 2006). We look therefore for solu- 
tions of ( [TtI i with uxb specified in the sense of (O. Using 
the ansatz A - A exp(/if -i-ife) with a generally complex A and 
any real k satisfying A; <K A;f we find first 

mik^ + s ft 

^ + (77T-e)A2jUJ" 



' A + {ri-Y + e)k^ 

m\k^ 



(27) 



The requirement of non-vanishing A poses an eigenvalue 
problem for A. The two eigenvalues, normalized to rj-ik^, are 



k^ 



1± — J(5A2+ 7712)7721 +e2. 



(28) 



A necessary and sufficient condition for an exponentially 
growing solution is that the radicand in ( |28] ) is positive and 
that it exceeds 77^. If the S term dominates and the others are 



neglected this condition turns into 

Sm\ 



(77tA)2 



> 1. 



(29) 



As k can be made arbitrarily small (by making the domain 
size large enough) this condition is always satisfiable if only 
57721 > 0. The neglect of the terms without S in the radicand 
is justified if I77127721 -I- e^l/77^ <k which can again always 
be guaranteed by sufficiently small k. Under this condition the 
maximum growth rate with respect to k is 5 772 1/4777 and oc- 
curs at A: = VS' 7721/2777. Consequently, as long as 7721 can be 
considered linear in S the maximum growth rate is propor- 
tional to and the corresponding k proportional to S . 

2.5. Simulations 

For the numerical simulations we use the Pencil Code*, 
where the test field algorithm has akeady been implemented. 
We employ periodic boundary conditions in the y and z direc- 
tions and shearing-periodic boundary conditions in the x di- 
rection (Wisdom & Tremaine 1988, Hawley et al. 1995) and 
use a resolution of up to 256^ meshpoints for the runs with 
the largest Reynolds numbers. As mentioned above, a com- 
putational domain of size is used, so the smallest finite 
wavenumber is Ai = Itt/L. As initial conditions for the hy- 
drodynamic part we assume vanishing velocity, U - 0, and 
uniform density equal to some value po- The initial condition 
in the test field calculations is a* - 0. Owing to the use of 
periodic boundary conditions, the total mass in the computa- 
tional domain is conserved, and therefore the mean density 
will be always equal to the initial value, (p) = po, where (. . .) 
denotes a volume average. 

In all investigations reported in this paper only weakly com- 
pressible turbulence has been considered. The Mach number 
Mrms/cs did not exceed a value of the order of 0.1. 

3. RESULTS FOR THE DIFFUSIVITY COEFFICIENTS 

As explained above, the test field procedure yields the co- 
efficients Tjij first as functions of z and t. However, after 
averaging over sufficiently long time intervals, we expect 
to approach the results for homogeneous, isotropic, mirror- 
symmetric and statistically steady background turbulence, in 
particular coefficients 77,^, being independent of z and f. We 
present here results for the 77,^ gained by averaging of the 
'raw' data first over z and then over time. In this context 
the effect of averaging over z consists in a first reduction of 

* |http : //www ■ nordita . org/sof tware/pencil- code| 
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the temporal fluctuations. This appears plausible in the pic- 
ture in which the domain contains a finite number of turbulent 
"eddies" (Hoyng 1993). We may interpret them as different 
realizations of a specific eddy and thus the average over the x, 
y and z of a given domain as an average over the ensemble of 
these realizations. When accepting the principle that the en- 
semble average is equivalent to a time average we see that the 
effect of averaging the original rjij over z is just equivalent to 
some temporal smoothing. After having averaged overz, time 
averages are then taken over a suitable stretch of the full time 
series where these averages are approximately steady. We use 
the time series further to calculate error bars as the maximum 
departure between these averages and the averages obtained 
from one of three equally long subsections of the full time 
series. 

Important control parameters that are being varied include 
the hydrodynamic and magnetic Reynolds numbers. Re and 
ReM, as well as the magnetic Prandtl number Pr^, with 

Re = Mrms/(v^f), RCm = Mrms/('7fcf), PrM = v/?? . (30) 

In the case of rotation we define further the Coriolis number 
Co and in the case of shear we define the parameter Sh, 

Co = 2Q/(M™s^f), Sh = Sl(u,^,kf) . (31) 

We note that Co, like Q, is never negative. In all cases with 
shear presented below, S and thus Sh are negative. For most 
of the calculations we use kfjki = 5, except in ^4. II where 
kflki - 10. In both cases the range of forcing wavenumbers 
is kf ± kill. 

3.1. Effect of rotation 

In the case of rotation (Co ^ 0), but without shear (Sh=0), 
the coefficients % and S are relevant. Figure [T] shows their de- 
pendence on the Coriolis number Co for fixed Reynolds num- 
bers. Re = 1.3 and Rcm = 13. We see that rjt shows a drastic 
decline when Co approaches and exceeds unity. This can be 
understood as a consequence of an evolving Taylor-Proudman 
state of the turbulent flow. Clearly S is positive. The ratio 6/ tji 
first increases with Co, but it begins to decline when Co has 
exceeded a value of about 3. In Fig. |2]the dependence of rj^ 
and S on Re^ is given for Re = 16 and Co = 1.3. As Re^ is 
increased, rj^ and 6 increase for Rcm < 10. 

According to the considerations in ^2.1l we have to expect 
that the diagonal elements 7711 and 7722 of the magnetic diffu- 
sivity tensor coincide. Indeed, the observed values of e (not 
shown) are only of the order of the errors. 

Our results are consistent with those obtained in the frame- 
work of the second-order correlation approximation, see Ap- 
pendix |B] We take this consistency as a confirmation of the 
correctness of the test-field method. 

3.2. Effect of shear 

We now discuss the case of shear (Sh 0) in the absence 
of rotation (Co - 0). Figure [3] demonstrates that the value of 
'7t/'7t{) clearly exceeds unity for not too small values of |Sh|, 
that is, shear leads to a slight enhancement of the turbulent 
magnetic diffusivity. At the same time, for negative values of 
Sh, both 7712 and 7/21 attain finite positive values. In Figs |4] 
and |5] these quantities are shown as functions of Re^, with 
Re = 1.4 and Sh - -0.6, or Pr^ = 20, respectively. 

We recall the dynamo condition ( |29] l. Since in all our simu- 
lations S is negative, a dynamo would be possible for negative 
7721 only. In Fig. |5] with Pr^ = 20, we see indeed negative 
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Fig. 1 . — Dependence of the normalized r]i and 5 on Co for homogeneous 
turbulence with rotation for Re = 1.3 and Re^ = 13. The vertical Hues 
indicate error bars. Note that there is a maximum of Sjrij at Co x 3. For 
Co < 3 the results for 6lrn are best described by Sjrn x 0.05 X Co""^', given 
by the dash-dotted line, but also a linear dependence, S/t]j x 0.035 X Co, 
indicated by the dotted line, is compatible within error bars. 
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Fig. 2. — Dependence of the normalized % and S on Re^f for homogeneous 
turbulence with rotation for Re = 1.3 and Co = 1.3. The dotted lines show 
the power law fits rjjrjto = 0.45 ReJJ^' and iS/j/t = 0.025 Re^ which apply for 
Rbm < 7. 



7721 for high Rcm- Considering the large error bars, however, 
we may hardly conclude that a dynamo is really possible. In 
general the errors could be reduced by extending the time se- 
ries. However, for large Re^ small-scale dynamo action oc- 
curs that introduces additional fluctuations whose amplitude 
increases exponentially with time, and so we have to stop the 
calculation. One remedy might be to reset f in regular time 
intervals, but this has not been done yet. 

As the considerations of ^2. II show, there is no general rea- 
son for an equality of the two diagonal elements 7711 and 7722 
of the diffusivity tensor, that is, e does not need to be equal 
to zero. As shown in Fig.|6] for Rcm of order 10 and above, 
e may deviate from zero, but its value is of the order of the 
error Again, our numerical results are in agreement with re- 
sults obtained in the second-order correlation approximation, 
see Appendix IB] 
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-Sh 

Fig. 3. — Dependence of T]t (normalized by i]to) as well as 7712 and (721 
(both normalized by rjj) for homogeneous turbulence with shear on Sh for 
Re = 1.4 and Re^ = 14. The dotted lines represent linear dependencies on 
Sh. 
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Fig. 4. — Dependences of rjt (normalized by ;;to), as well as rji2 and (721 
(both normalized by t]j) for homogeneous turbulence with shear on Re^ for 
Re = 1.4 and Sh = -0.6. The dotted lines con'espond to rjjriio = O.SReJJ^*, 
Vll/VT = 0.2 InRcM and (721 /'?T = 0.05 InRe^ and illustrate that rmlryi and 
'721 /J^T vary only weakly with Re^. 



COMPARISON WITH DIRECT SIMULATIONS AND SIMPLDTED 
MODELS 




Fig. 5. — Same as Fig.|4] but for a fixed magnetic Prandtl number, PrM = 20. 
Here, Re = ReM/20 is not constant. Because u^^^ increases with increasing 
Rbm, Sh is also not constant and varies between -2.5 (for Re = 0.16) and 
-0.3 (for Re = 13). Note that 7721 turns negative at about Rcm = 100. How- 
ever, the errors are larger than the mean. 
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Fig. 6. — Dependences of e (normalized hyifi) on Sh (left) and on Re^ 
(right) for the same runs as in Fig.[3]and Fig.|4] respectively. 



4.1. Large-scale fields in simulations with shear 

We report now on calculations with the original induction 
equation ( fTSl l, instead of the test field equations ( l20l l. together 
with the hydrodynamic equations (|9]l and ( fTOl i. In the momen- 
tum equation however, the Lorentz force was restored, 
thus providing a nonlinear feedback of the magnetic field. In 
all cases we used Sh = -0.15 and kflk\ - 10. The results 
are shown in Figs [T] and [8] for two different combinations of 
ReM and Pr^- In both cases there is an initial phase where 
the mean field grows exponentially. Mean fields with a par- 
ticularly prominent B, component occur. The B ^ component 
seems to be in antiphase with By, as expected for negative 
shear, but this component is much more noisy. Furthermore, 
for Pr^ = 7 and Re^v/ - 130 (Fig.|7J there are episodes where 
By fades away and later reemerges, but possibly with the op- 
posite orientation. Similar results (not shown here) have also 
been obtained for smaller values of Rcm- In the case with 
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Fig. 7. — Time dependence of the rms value (with respect to z) of By (top) 
and space-time diagrams Bv(z. ') and By(z, t) (all in units of Beq. where Beq = 

Vpo(pw^)) from a direct simulation with Re^ = 130, Pr^ = 7, kijki = 10 
and Sh x -0.15. The top panel demonstrates the initial exponential growth 
of the mean field (the growth rate is 0.009 Mrms^f)- The other panels show 
episodes of large scales in z especially in the By component. 

PiM - 20 and Re^ = 200 (Fig. [SJ, however, B,, keeps the 
same orientation throughout the run. 

4.2. Magnitude and effect of fluctuations 

In ^3.21 we have seen that the sign of 7721 is not suitable 
for enabling a shear-current dynamo except perhaps for high 
values of Rsm- On the other hand, as demonstrated in ^4.11 
large-scale magnetic fields are being generated. Explaining 
this in terms of the shear-current effect is very questionable. 
Therefore we ask now whether an incoherent alpha-shear dy- 
namo (Vishniac & Brandenburg 1997) might play a role. An- 
other explanation would be an incoherent shear-current dy- 
namo that we discuss below in ^4.31 

The possibility that a random a with zero mean can pro- 
duce magnetic fields was first discussed by Kraichnan (1976) 
and Moff'att (1978). In the presence of shear, strong large- 
scale fields can be generated (Vishniac & Brandenburg 1997, 
Sokolov 1997, Silant'ev 2000, Fedotov et al. 2006, Proctor 
2007). Consider an incoherent alpha-shear dynamo with a 
scalar a fluctuating around zero. In the limit A; <sc fcf and if 
\a\k «: \S I the condition for mean fields growing on the aver- 
age exponentially reads 

£>.S =t^rms|5 1/(772^3) >DS', (32) 

where D"" ~ 2.3 for a white noise a eff'ect; see Appendix ICl 
In a finite domain all mean-field coeflicients show fluctua- 
tions, and so a/y must fluctuate about zero. We may extend the 
test-field procedure for the determination of the 77,^ such that 
it provides us the atj, too. When starting from dU instead of 
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Fig. 8. — Same as Fig. [7] but for Rsm = 200 and Pr^ = 20. Initially, the 
field grows exponentially at a rate 0.012 Mnyn^A'f. Note that the mean field is 
nearly steady. 
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Fig. 9. — Dependences of the rms values of the temporal fluctuations am,, 
(normalized by ;;t0^f), 'y"™' ^'\2 (normalized by );to), on Re^ for 
Re = 1.4 and Sh = -0.6. 

(O and using again the four test fields B^'^ and B'*^, q - 1,2, 
we find 

an = B-^{&Y cos kz + I-Y si" fe) 

aa^B-^{Efco^,kz + 8f?.mkz), /=1,2, (33) 

together with the relations ( |23] ) for 77,1 and 77,2 (see also Bran- 
denburg 2005b). 

In contrast to the considerations in Owe consider now the 
mean-field coeflicients, as obtained from the test field calcu- 
lations, after averaging over z, but not over t. Then the aij 
consist of fluctuations around a zero mean, that is, we have 
an incoherent a effect. (Without the averaging over z the fluc- 
tuations would be even bigger) In the case of fluctuations of 
7721 and ?7i2 we speak analogously about an incoherent shear- 
current effect. 
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Fig. 10. — Probability density functions (PDFs) of Q;,y (top) and 77,^ (bot- 
tom) for Rbm = 14, Re = 1.4, and Sh = -0.6. The PDFs for the different 
components of a/y (dotted lines) are close together; their average is given by 
the solid staircase line and compared with a Gaussian fit. The PDF of rj2\ 
(solid line) is around zero while those of ;;i2, f/u, and rjii are not. (The latter 
two are simply denoted by rjt.) 

We have calculated the rms values of the temporal fluctu- 
ations of the ajj and 77 ,y which are denoted by a™^ and 77™\ 
respectively. We have also taken the averages over all four 
components of a™^ and over the two diagonal components 
of 77™'* and denoted them by o-™^ and 77™% respectively. Fig- 
ure |9] shows these quantities along with 77™^ and 77™** for a 
small Reynolds number and moderate shear as functions of 
Rcm- They are all mildly growing. 

In Fig.[TO]we show that the probability density functions of 
aij and 77,^ for a run with Re^ = 14, Re = 1.4 and Sh = -0.6 
are approximately Gaussian. In order to improve the statistics 
we have, in addition, averaged the results for all four compo- 
nents of ajj. The result is similar to those for the individual 
components. The diagonal components of 77,7 are distributed 
around finite averages, and 7721 is distributed around a positive, 
but small value. 

4.3. An incoherent shear-current dynamo? 

Yousef et al. (2007) have reported large-scale dynamo ac- 
tion at low Reynolds numbers (Re = Re^ - 5) for weak shear 
[Sh < (37r)"'] in tall boxes so that the smallest wavenumber 
in the z direction, k\-, can be up to 128 times smaller than 
those in the other two directions. They discuss in more detail 
the case where it is 16 times smaller, i.e. ki- - ki/16. Using 
small Reynolds numbers has the advantage that small-scale 
dynamo action is then impossible. 

We have analyzed similar cases (Re - 5 with Pr^ = 1 and 
0.01 < -Sh < 0.3) using however cubic domains of size L^, 
so kiz = ki, and a forcing with kf/ki equal to 5 instead of 3. 
It turns out that the value of the crucial coefficient 7721 fluc- 
tuates around zero. This is also plausible from Fig. [3] (even 
if it does not apply to Ptm = 1). We must therefore con- 
clude that the (coherent) shear-current eff'ect cannot explain 
the generation of the mean magnetic field found by Yousef et 
al. (2007). With respect to the incoherent effects it can be seen 
from Fig.[TT]that the values of a™'^ and 7721"* are more or less 
the same for different Pr^, Re^, and Sh. 

Let us consider the instantaneous values of the growth rate 
A, as calculated from equation ( l28l l with the fluctuating 77,^. If 
k = ki, we get always negative A. However, for the smaller 
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Fig. 11. — Dependence of the rms values of the temporal fluctuations ff'"" 
(normalized by r/to^f) and r/'J^^ (normalized by rjto), on Sh for models with 
PiM = 1 and Re = 5 (dashed lines) compared with the models shown in 
Fig.[3]with PrM = 0.1 and Rsm = 14 (solid lines). 

value k - ki/16, appropriate for the model of Yousef et al. 
(2007), it is possible to have large positive A during extended 
periods of time. Although B can be amplified during those 
episodes it must decay during episodes with the reversed sign 
of 7721, and it is not certain from this competition whether a 
dynamo powered by the incoherent shear-current effect may 
result. 

So far we ignored the possibility of an incoherent alpha- 
shear dynamo that must work at the same time. In order to 
assess the relative importance of the two incoherent effects we 
have considered a simple model with random a and 77 tensors 
that are delta-correlated in time. (Delta-correlated noise is 
the simplest model; a more realistic case would be to assume 
colored noise with a finite correlation time.) The model is 
explained in Appendix ICl Under the assumptions \ai\\ <K 
15 1 /A: and |q'i2|, \a2i\ ^ '/t^, as well as I7712I <s; \S\/k^ and |e| <K 
|77t|, its governing parameters are the two dynamo numbers for 
the incoherent effects. 
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which can be expressed in terms of the three quantities 
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In Fig. [T2]we give a contour plot of the normalized growth rate 
as a function of the two dynamo numbers, Das and D^is . For 
small values of , the incoherent shear-current effect has 
a slightly adverse effect on the dynamo, but for larger values 
it lowers the critical value of Das significantly. For Das - 
even a purely incoherent shear-current dynamo is possible if 
D„s > 6.5. 

In interpreting simulations, we focus on domains whose 
smallest finite wavenumber in the z direction is kir. In view of 
the already discussed dynamo found by Yousef et al. (2007) 
we employ the data for a™'* and 77™^ originating from our 
aforementioned calculations of similar cases to derive growth 
rates with the help of Fig. [12] With k - k\^ - ki/16 we find 
first Cs ~ 40, and from Figs |5] and [TT] we have Q ~ 0.1 
and Ca ~ 0.1 (using kf/ki = 3), so that D,js ~ 4 and 
D„s ~ 4. This suggests that the incoherent shear-current dy- 
namo is subcritical, while the incoherent alpha-shear dynamo 
is supercritical. Therefore an incoherent alpha-shear dynamo 
seems to be a plausible explanation. This explanation is fur- 
ther supported by our finding that for constant rms values the 
growth rate is in good approximation a linear function of 5 - 
just as observed by Yousef and coworkers. 
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Fig. 12. — Contour plot showing the normalized growth rate as a function 
of the dynamo numbers for the incoherent alpha-shear and shear-c urren t dy- 
namos from a numerical solution of the model described in Sect. 14.31 The 
zero line is given in white. 

An explanation in terms of the incoherent a effect is also 
suited for the nonhelical dynamo of Brandenburg (2005a), 
where Cs ~ 25, C,, ~ 0.15 and C„ ~ 0.2 for the appropri- 
ate value of the magnetic Reynolds number, Rsm = 80, so 
that Dijs ~ 4 and Ars ~ 5. For this model we used kf/ki = 5, 
although the shear profile is here more complicated. 

Finally, let us return to the cases of dynamo action consid- 
ered in Figs|7]and|8] If we assume that the negative values of 
7721 seen in Fig. |5]for Rsm > 100 are real, we have to ask for 
the relative importance of the regular (coherent) shear-current 
effect and the two incoherent effects. For both dynamo cases 
in Figs I2] and |8] the values of Cs are ~ 30 (but uncertain be- 
cause the dynamo-generated field quenches the value of rjt), 
so the three dynamo numbers would be D,js ~ 1-5 for the 
regular shear-current effect, D,js - 1 .5 ... 3 for the incoher- 
ent one and D^s = 3 ... 6, respectively, where we have used 
= 0.05... 0.1 and C„ = 0.1 ...0.2 (cf. Figs|5]|9]and[n]l- 
The values of D„s and could be somewhat smaller if one 
takes into account that the level of fluctuations is smaller for 
kf/kiz - 10 instead of 5. We recall that the corresponding 
critical values are 1, 6.5 and 2.3, respectively. Hence, with 
respect to the regular shear-current dynamo this case is only 
slightly supercritical, but subcritical with respect to the in- 
coherent shear-current effect and supercritical with respect to 
the incoherent alpha-shear dynamo. By inspection of the val- 
ues of the growth rate it is possible to infer safely that this 
situation is dominated by the incoherent alpha effect. The in- 
coherent shear-current effect has a weakly adverse influence 



whereas its regular counterpart clearly supports dynamo ac- 
tion. 



5. CONCLUSIONS 

The present work has demonstrated that the test field 
method provides a robust means of determining all compo- 
nents of the turbulent magnetic diffusivity tensor that are rel- 
evant for mean fields depending only on z and t. Both ro- 
tating and weakly shearing turbulence are studied. In either 
case the diagonal components of the turbulent diffusivity ten- 
sor are about equal to each other. Shear slightly enhances the 
turbulent magnetic diffusivity while rotation quenches it. In 
the presence of rotation, the ilx J effect occurs, which is de- 
scribed by the off-diagonal components of the turbulent mag- 
netic diffusivity tensor Shear leads to the shear-current effect, 
again described by off-diagonal components of this tensor In 
both cases the results are consistent with those found in the 
framework of the second-order correlation approximation. 

The possibility of the so-called shear-current dynamo has 
been scrutinized. It depends crucially on the sign of the com- 
ponent 7721 of the magnetic diffusivity tensor. It turns out that, 
within the ranges of parameters considered, its sign is in gen- 
eral not suited for driving a dynamo based on this effect, with 
a possible exception at large magnetic Reynolds numbers. In 
this way the analytic results found in the second-order cor- 
relation approximation for incompressible fluids (Riidiger & 
Kitchatinov 2006; Radler & Stepanov 2006) are confirmed 
and generalized. 

Direct numerical simulations are presented which exhibit 
growing mean magnetic fields in shear flow turbulence. An 
interpretation as a (coherent) shear-current dynamo is hardly 
possible. Instead, it is argued that it can be explained by an in- 
coherent alpha-shear dynamo. The incoherent shear-current 
effect has also been determined, but it is found to be less im- 
portant. 
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APPENDIX 
CONCERNING = 

In the case of rotation the tensorial structure of q-^ must agree with that of rjij given in (O, that is, 

aij - ao6ij + aieijkClk + a20.iQ.j . (Al) 

Since aij is a pseudo-tensor and Q an axial vector, the coefficients ao, ai and q'2 must be pseudo-scalars. Under our assumptions, 
however, no pseudo-scalars can be constructed. So we have to conclude that aij - 0. 
In the case of shear we may argue analogously. Referring to (|5]l and the subsequent explanations we have then 

aij = aodij + aigigj + a2hihj + a^gihj + a/^gjhi (A2) 

with ao, a\ ■ ■ ■ a^ being pseudo-scalars. Again, it is impossible to construct pseudo-scalars. Thus, we have again a,/ = 0. Of 
course, the situation would be different if the shear provided (large-scale) kinetic helicity, as then the pseudo-scalar ■ curl 
would be available. 
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Fig. 13. — Results for an incoherent alpha-shear dynamo (D,,^ = 0). Left: Temporal evolution curves of Bnns showing that the critical value of D^s is around 
2.3. Nonlinearity is here ignored and Bin,s is scaled by its initial value, Bo. Right: With additional dynamical quenching, D^s = 50 and k[/ki = 5. Note that the 
typical time between reversals increases with Re^ approximately like ReJ^^. 



COMPARISON WITH RESULTS OF THE SECOND-ORDER CORRELATION APPROXIMATION 

In a paper by Radler & Stepanov (2006, referred to as RS06 in the following) the mean electromotive force has been calculated 
in the second-order correlation approximation for generally inhomogeneous turbulence in an incompressible rotating fluid show- 
ing a position-dependent mean motion. In this context the second-order correlation approximation was understood as the neglect 
of higher-order terms in the induction equation as well as in the momentum balance. Both the Coriolis force and derivatives of 
the mean velocity were assumed to be small enough so that the mean electromotive force is linear in the angular velocity Q and 
the gradient tensor of U. Detailed results were obtained for a special correlation function of the background turbulence. 

Let us apply the results to the situations considered in the present paper. In the case of rotation without shear we obtain 

- = 7 McoRsM(A,kff ^/^6^q,PrM) (Bl) 
7]t 4 V 2 

with Co, Rqm and Pr^ as defined above, q = A^jrjT^, and /Ic and Tc being correlation length and time, respectively. When 
introducing the Strouhal number St = Mrms^f^c^ we have q - Re^/St. It seems plausible to assume that A^kf » 2n. The function 
(5° can be calculated according to 5° = [5"*"^^, Prw) + /^°*"^(^, PrM)]/2^'(°)(^) from the functions (5°(">, and defined 
and plotted in RS06. It turns out that 6^ is never negative and approaches unity if Pr^ = 1 and ^ — » 0. Of course, we have 
^hi - {5 hdixithi)- The factor rjJrj-Y depends on Re^, /ic^f and 0^^^Hq). It satisfies < tjxItji < 1 and approaches unity as 
Rcm — » 

Clearly (IB 11 1 and the results reported in 33. II agree in the sign of 6. Although these results do not really confirm the linearity 
of 5 in Co, which is suggested by (IB 11 1, they are not in conflict with that, see Fig.[T] A further comparison of results is difficult 
because of, e.g., the not exactly known value of and the errors of the data presented above. 

Proceeding to the case of shear without rotation we note first that, due to the aforementioned assumption on the linearity in the 
mean- velocity gradient, that is in 5, both k\i and K22 are equal to zero. Furthermore, we have 

^ = -\shRtM{A,kffif,^{q,?yM) , — = ~ShRsM(A,kff4^(q,PrM) . (B2) 
?7t 5 Tjt 5 

Here = i[77"(^>(^,PrM) + 77"<"''(?,PrM)] and 77^; = i[77"(^'(^,PrM) - ;7"<**''(?,PrM)], where t;"*^* = [13A:''(^^(^,PrM) - 
7^0(«)(^,PrM)]/6yS»<»>(<7) and /yOf^' = [550<"'>(<?, PrM) + /^''<^>(^, PrM)]/6/?"("'(^), with the functions k'^'-'^K j^'^'^K 5°^, ^^^^ and 
^O(O) RS06. The quantities rf!^^ and rf^^ approach unity and zero, respectively, if Pr^ = 1 and ^ — > 0. We note that -rjnlS and 
7721/5 coincide with the quantities 6' and S introduced in Appendix D of RS06, respectively. It has been shown there that this 5 
(different from that considered above) cannot take negative values. This applies then to 7721, too. 

Being aware that the second-order approximation applies only for Re^ that is not too large, we may state that (IB2l i and the 
numerical results reported in 33.2l agree in the sign of 7721. The possible deviation in Fig.|5]is outside the validity range of this 
approximation. The linearity of 7712 and 7721 in Sh indicated in ( IB2b is well confirmed by the numerical results; see Fig. |5] Again, 
further comparison of the results is, for the reasons mentioned above, rather difficult but no striking disagreement has been found. 
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INCOHERENT ALPHA-SHEAR AND SHEAR-CURRENT DYNAMOS 

We calculate numerically solutions of the dynamo equation with incoherent alpha and shear-current effects in unbounded 
space. It reads 

DA _ _ _ 

— ^-SA,x + G-j]J, (CI) 

where _ _ _ 

8i = aij(t)Bj-T]ij(t)Jj. (C2) 

The a and r] tensors are delta-correlated in time with 

(aijit) a,/f')) = (<f 6(t - t'), {aij) = 0, (C3) 

(rjijit) riijit')} - {jjijf = (77™f 6(t - t'), (C4) 

where no summation over double-indices is assumed and (. . .) means here a temporal or ensemble average. For solving ( IClb we 

use the ansatz 

A(z,f) = A(f)exp(ife) (C5) 

with an arbitrary, but fixed wavenumber k and employ a third-order Runge-Kutta time stepping scheme. At each time step of 
length 5t, the fluctuations of aij and rjij are taken as random numbers from a Gaussian distribution and scaled by 1 / ^[6t so that 
equations ( |C3] ) and ( IC4I ) hold. 

We recall that, if aij — adij with a - const S /k, the critical value of the dynamo number Das as defined in equation ( |32] | 
but with Q'"™ replaced by a, is D™' = 2 (e.g. Brandenburg & Subramanian 2005a). In the case of a pure incoherent alpha-shear 
dynamo, i.e. D^s - 0, it is found that D™' ^ 2.3 (see the left panel of Fig. \T3[ . On the other hand, for D^s = we have 

There are reversals on a typical timescale of about one diffusion time. However, this time can increase significantly if magnetic 
helicity conservation (appropriate for a closed domain) is taken into account (Field & Blackman 2002, Blackman & Brandenburg 
2002, Subramanian 2002). This means that the a effect has to be amended by an additional term that results from the current 
helicity produced by the dynamo. We assume again aij - adij and a non-fluctuating rfij - ritSjj, further 

a(t) = aK{t) + auit), (C6) 
where q'k(0 is stochastic, just like a,/?) in equation iC3i . and Q'M(f) obeys the differential equation 



-2T]^k^ 



df 



Re 



M 



(C7) 



with B and fi defined analogously to A in (IC5b . The dynamo number D„s is now defined with respect to a™^ . Equation ( IC7I ) is 
solved simultaneously with equation dClb using the aforementioned time stepping scheme. As here a nonlinearity is introduced, 
the ansatz ( IC5I ) has now to be understood as a one-mode truncation. The model calculations show that the timescale for reversals 
increases proportional to Re (see the right panel of Fig. [T3] l. 
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